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Magnetic field reversals and galactic dynamos 
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We argue that global magnetic field reversals similar to those observed in the Milky Way occur quite frequently in mean-field galactic 
dynamo models that have relatively strong, random, seed magnetic fields that are localized in discrete regions. The number of reversals 
decreases to zero with reduction of the seed strength, efficiency of the galactic dynamo and size of the spots of the seed field. A systematic 
observational search for magnetic field reversals in a representative sample of spiral galaxies promises to give valuable information 
concerning seed magnetic fields and, in this way, to clarify the initial stages of galactic magnetic field evolution. 
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1 Introduction 

An outstanding question in the study of magnetic fields in spiral galaxies is that of large-scale field reversals. 
Do they occur, how easily can they be detected, can dynamo models explain them? 

It appears that there may be such a field reversal in the Milky Way - certainly there is a localized 
feature relatively close to the Sun, but the difficulties of observing through the confusion of the disc plane 

aimed 



1996 



make it difficult to be absolutely certain wheth er it is global in extent. Further reversals have been c 
to exist, but with increasing uncertainty fe.g. ISimard-Normandin and Kronberg 19791. Beck et al. 
Frick e_ alJboOll . iMen et a/Jl20o3 . iNota and Katgertll2010l . IPshirkov et alJl201ll . Ivan Eck et q/.ll2011 
201 lh . Speaking generally, the idea that magnetic field reversals are present in the Milky Way is widely 



Beckl 



accepted. 

The situation in external galaxies (see the review Beck et al. 1996), where the problem of observing from 



a position in the disc plane is absent, is nevertheless rather more uncertain. Except for the immediate 
neighbours of the Milky Way, the available resolution is too low to make definitive statements. But in the 
well observed nearby M31, reversals appear to be absent. 

Classical mean-fi eld dynamo models that start from a dynam ically small field do not in general produce 
field reversals (e.g. Ruzmaikin et all 19881 . Meinel et al\ 1990T ). because the leading eigenfunction of the 



kinematic mean-field galactic dynamo usually has no reversals, and this dominates the field structure as 
the field strength grows. 

The situation becomes less straightforward when the seed field for the dynamo is stronger in comparison 
to the contemporary galactic magnetic field, so the galactic dynamo becomes nonlinear before the magnetic 



these are essentially fiel d discontinuities, (see e.g. 



(Poezd et ollll993l.lBeck et al 1994. 


Moss et al 1998 V. 


Belvanin et al. 


1994. 


Vasil'eva et al. 1994. Moss et al. 



existing examples are one dimension al. 
More recently, Moss et al. ( 20121 ) have published an extended mean-field model in which long-lived 



reversals appear to be a quite common feature. This dynamo model differs in several important ways from 
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earlier 2D mean-field models. Notably these computations begin from initial conditions of small-scale, 
approximately equipartition strength, fields that are randomly distributed in many discrete spots. This is 
in contrast to nearly all previous mean- field models, which begin from dynamically weak seed fields, either 
quasi-homogeneous or small-scale. The rationale is that these dynamically strong fields are the result of 
small-scale dynam o action in star forming regions in the protogalaxy. In the computations described in 
Moss et al. ( 20121 ) there are also ongoing injections of small-scale field, but we will argue here that these 
are not central to the generation of reversals. 

Our approach is close to that of the one dimensional model of lBeck et a/1 (|l994l ) , where a localized random 
initial condition was considered and long-lived large-scale magnetic reversals developed. (This contrasts 
with the model oflPoezd et al. ( 1993), wh ere there were ongoing injections of random fields.) Thus both 
Poezd et al. (|1993T ) and iBeck et al\ ll994h considered magnetic fields depending on galactocentric radius 
only, i.e. ring-like magnetic structures were described. The novelty of our approach is that we consider seed 
fields that depend on both radius and azimuth, i.e. there is no a priori restriction to ring-like structures, 
and any such structures develop as a result of the dy namo proces s . We demonstrate below that ring-like 
structures with reversals, similar to those described b vlBeck et al. (1994). also develop in our simulations. 



Thus our model reinforces and extends the results of Beck et al 



"(jl994T ). 



On the other ha nd, we note that the stable magnetic ring-like structures with reversals obtained by 
Poezd et al\ (j 19931 ) appeared in a model with some fine tuning of the model parameters. Stability conditions 
for such structures obtained bv lBelvanin et al\ (j 19941 ) depend on rather subtle properties of the galactic 
disc; this gave a hint that the magnetic configuration of the Milky Way with field reversals is a rare 
exception to the typical situation without reversals. Possibly, there was a misinterpretation of the results, 
and maybe the authors of the papers cited above were not insistent enough. Nevertheless, dynamo generated 
configurations with reversals remained rather neglected by the astronomical comm unity. Moreover, th e 



option of including continuous injection of small-scale magnetic field as considered by Poezd et al. (jl993l ). 
was not further developed until recently. 



In contrast, the new generation of dynamo models presented in iMoss et al\ (120121 ) produce magnetic 
structures with reversals even for the simplest and most primitive distribution of the dynamo governing 
parameters and they appear to occur more or less as commonly as structures without reversals. This is 
why we consider below examples of dynamo models with (and without) reversals. They are demonstrably 
oversimplifie d and do no t inclu de any fine tuning of the governing parameters of the dynamo. 

The paper IMoss et all (12012 ) focussed on applications of the model to future observations of magnetic 
fields of the earliest galaxies using the recently developed (e.g. LOFAR) and planned (SKA) radio tele- 
scopes. The aim of this paper is to describe the long-lived reversals obtained as a quite general phenomenon 
of galactic dynamos. We are motivated here by future observations of magnetic fields in very young galax- 
ies, whose detailed structures are almost unknown a t the momen t, and so we study a very sim ple, generic 
model . This is in some contrast to the approch of Poezd et al. ( 19931 ) (and to some extent IBeck et al. 



( 19941 )) who focussed their attention on the dynamo properties of particular nearby galaxies (M31 and the 
Milky Way) and exploited particular subtle properties of their rotation curves and other parameters. 



2 The dynamo model 



Detailed modeling of galactic magnetic of galactic magnetic field evolution from the time of galactic 
formation up to the age of a contemporary galaxy is a severe problem, both be cause of oby i ous nu merical 



difficulties and our very poor knowledge of hydrodynamics of early galaxies. IMoss et al 



t()l$ ) used a 

reasonable simplification of the problem in the mean-field approach in the form of the "no-z" model (e.g. 
Mosd 1995 ) which restricts the modelling to quantities which are accessi ble obseryation ally, at least in 
the immediate future, and make the numerical implementation afford able. Phillips! d2001f) sy stematically 
compared this approximation with the standard "local model" (cf. kuzmaikir, et all Il988h and found 
satisfactory agreement. 

He re we use the same simplified model. For the sake of clarity we briefly reproduce the relevant equations 
from IMoss et all (12012 ). The code solves in the auo approximation explicitly for the field components 
parallel to the disc plane while the component perpendicular to this plane (i.e. in the z-direction) is given 
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by the solenoidality condition. An even (quadrupole-like) magnetic field parity with respect to the disc 
plane is assumed. The field components parallel to the plane can be considered as mid-plane values, or 
as a form of vertical average through the disc. The key parameters are the aspect ratio A = h/R, where 
h corresponds to the semi-thickness of the warm gas disc and R is its radius, and the dynamo numbers 
Ra = ctoh/T],R u = f2oh 2 /r] (if disc thickness varies with radius then h is some reference value, ho say). A 
must be a small parameter. r\ is the turbulent diffusivity, assumed uniform, and ao, Oo are typical values 
of the a-coefficient and angular velocity respectively. Thus the dynamo equations become in cylindrical 
polar coordinates (r, (f>, z) 



~gf = ~R a B^ - —B r + A 2 



9 f 1 9 ( R 

or \r or 



+ 



1 d 2 B r 2 dB d 
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1 d 2 B d 



2 dB r 



(2) 



where z does not appear explicitly. Here the disc semi-thickness h is assumed to be constant: taking 
h = h(r) would introduce geometric correction factors (h(r) / ho) , (h(r) / ho) 2 multiplying R a , and the 
terms — \{B T , B^), respectively. This equation has been calibrated by introduction of the factors 7r 2 /4 in 
the vertical diffusion terms. In principle in the acv approximation the parameters R a , R u can be combined 
into a single dynamo number D = RaR^, but we choose to keep them separate. Equations (1) and (2) are 
displayed in polar coordinates but for ease of numerical implementation they are rewritten in Cartesians. 
Computations are performed on a unifo r m gri d of resolution 229 x 229 in the region —1 < x,y < 1. 
This resolution was tested in iMoss et all (|2012h . and indeed is in excess of what is required for the very 
straight forward dynamo p roblem solved after t = 0. The unit of time is about 0.8 Gyr. Further details are 



given in Moss et al. ( 20121 ). 



Length, time and magnetic field are non-dimensionalized in units of R, h 2 /rj and the equipartition field 
strength B eq respectively. A naive algebraic a-quenching nonlinearity is assumed, a = ao/(l + B 2 /B 2 q ), 



'eq 

where B eq is the strength of the equipartition field in the general disc environment. Note that small-scale 
field is injected only at time zero; after that we conduct a standard mean field dynamo simulation. 

The coefficients ao and B eq are assumed to be uniform in the work described in this paper because 
our restricted knowledge of the properties of the earliest galaxies give no secure basis for more realistic 
assumptions. For example, it is argued that «o depends on both disc thickness and local angular velocity, 
i.e. ao oc Qh~ l , where h can be expected to increase with radius, and Q to decrease. Of course, we continue 
our simulations to a time corresponding to the present day where knowledge of particular galaxies is rather 
better - we emphasize again that we are studying generic properties of thin disc dynamos. Taking typical 
galactic values, we can estimate R a = 0(1), R^ = O(10) (i.e. D = O(10)). 

We appreciate that more sophisticated approaches to the saturation process exist, however it does appear 
that in some cases at least, such a naive algebraic quenching can reproduce reasonably well the results 
from a more sophisticated treatment (e.g. iKleeorin et aZ.ll2002l ) . 



3 Results 



Moss et a/] (l2012h presented a particular galactic dynamo model that quite often resulted in magnetic 
field config u ration with one or more reversals. We argue here that many specific features of the model 
Moss et al\ (|2012h are not essential for the development of reversals. 

We show in figure the s tatistically steady field (at age ca. 13 Gyr) resulting from a standard com- 
putation of IMoss et aU (|2012l ) (model 135, R a — 1, R^j — 20; we keep here and below the numbers of the 
models from the working journal to allow identification of the models between various papers). In this case 
the initial field is randomly distributed in n spo t = 100 discrete regions, with r.m.s. value -Bo of the field 
in these regions of order unity - i.e. the alpha-quenching nonlinearity is immediately important. Figure Q] 
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Figure 1. (Colour online). A typical configuration of the initial magnetic fiel d; only the B vect ors at the centre of the spots are shown. 
The shaded circles indicate the gaussian half-width of the spots. (Taken from Moss et al\ j2012l U 



shows a representation of this initial field. Large-scale reversals are clearly visible in the resulting statis- 
tically steady configuration of this model - figure [2^,. We then show in figured (model NOINJECT10) 
the results of a similar computation using the same initial conditions, the only difference being that there 
is now no ongoing injection of fields in discrete regions. The final, steady, field is now smooth (the irreg- 
ularities visible in figured are caused by the ongoing injection of small-scale field), but two large-scale 
reversals are again clearly visible. A reversal has appeared by dimensionless time t = 1 (age < 1 Gyr), 
and have stabilized by t = 3 (ca 2.5 Gyr). A time sequence is shown in figure[3]for the model whose steady 
configuration is shown in figure [2)3. 

We expect that by reducing the strength of the seed field we will reduce the number of reversals to zero. 
In order to test this prediction we repeated the computation of figure [2)3 with the r.m.s. value, Bq, of the 
initial field in the discrete spots that was smaller by a factor of 10 -6 , and then by a factor of 10~ 8 . In the 
first of these cases the steady configuration has a single reversal (figure U^, model NOINJECT11), in the 
second case there are no reversals (figure Ht, model NOINJECT16). 

We also anticipate that a reduction of efficiency of the large-scale dynamo will affect reversals in a more 
or less similar way as reducing the seed field strength. The point is that both of these changes act to 
prolong the stage of quasi-kinematic magnetic field evolution. We verified this expectation, repeating the 
above computations with R u reduced from 20 to 10; then with Bq ~ 1 there is one reversal in the final 
steady configuration, and with R u = 7 there are no reversals. With B « 10~ 6 , R u = 20 there are none. 

If the size of the spots that localize the seed field becomes very small, reversals disappear. We demonstrate 
this in an experiment where the initial field is randomly assigned at every mesh point, again with R a = 
1, R u = 20, and Bq ~ 1. The final steady field configuration is shown in figureEb (model NOINJECT12) - 
it is smooth with no reversals. A similar result was found with Bq ~ 10~ 6 . Note that, in order to simplify 
comparison between these models, we used the same realization of random numbers to determine the seed 
field distribution in each case discussed above, with the exception of the case shown in figure SJa. 

We further show that the eventual steady state configuration has some dependence on the initial condi- 
tions, by the follo wing experiment. The case NOINJECT10 (the same initial conditions as model 135 from 
Moss et al. (120121 ). without the ongoing field injections) is re-simulated with a different realization of the 
random numbers, arriving at a field configuration that is quite different from that shown in figure [2}} - see 
figure 0Ji. We also repeated the NOINJECT10 calculation with 'global' alpha-quenching, i.e. quenching by 
the mean global energy. The final configuration was quite different, but reversals are again present. We 
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Figure 2. (a) Statistically steady state of model 135 of Moss et al. (2012) with R a = 1,Roj = 20 at dimensionless time t = 17 
(approximately 13 Gyr); (b) steady state of similar model without the ongoing injections of random field (NOINJECT10) at the same 
time. 



(a) 
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(b) 



(c) 




Figure 3. Snapshots of the field configuration of the model NOINJECT10 at times (a) t = 0.1 (ca. 0.08 Gyr); (b) t = 1.0 (ca. 0.8 Gyr); 
(c) t = 3.0 (ca. 2.4 Gyr). 

conclude that the form of the initial field can influence significantly the steady state configuration. 

Repeating the case illustrated in figure 2b with n spo t = 50 produces a superficially quite similar final 
state. As a final check, a simulation with n spot = 1, Bq = 1 produces a field without reversals. 



4 Discussion and Conclusions 



We deduce from these experiments that long-lived reversals can be generated by a standard 2D (nonax- 
isymmetric) mean-field dynamo model, given suitable initial conditions and a suitable range of the dynamo 
governing parameters. The initial conditions play a key role - this is probably why such solutions have not 
been widely noted previously. Different steady solutions can be found for the same dynamo parameters, 
depending on the initial conditions. Large-scale reversals can appear after about 1 Gyr, and they stabilize 
after 2 — 3 Gyr. The important elements appear to be an inhomogeneous distribution of the initial field 
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Figure 4. Steady field configurations for (a) R a = 1, R^ = 20, r.m.s. value of initial spotty field Bo ~ 10~ 6 ; (b) R a = 1, R^ = 20, initial 
field randomly assigned at every mesh point with r.m.s. value of unity; (c) as (a) with Bq £S 10 — 8 ; (d) the model NOINJECT10 with a 
different realization of random numbers - compare figure [2p. 



and strong differential rotation (larger Ru), the latter effect being enhanced by initial fields that are locally 
near equipartition values. Note that even with larger values of (R^ = 20 say), some initial conditions 
still produce final configurations without reversals. 

We note that some of these results were anticipated by Beck et all 1 19941 ). This paper studied a one- 
dimensional model and showed that, with a sufficiently strong small-scale seed field, reversals could persist 
for times of order the age of the Universe. 

As a straightforward consequence of the above conclusion, a systematic observational search for magnetic 
field reversals in a representative sample of spiral galaxies may provide valuable information concerning 
seed magnetic fields, and so clarify initial stages of galactic magnetic field evolution. This appears to be 
a realistic possibility for the forthcoming generation of radio te lescopes. For a contem porary comparison, 
it might be appropriate for orientation to look at figure 11 of van Eck et al\ (1201 ll ). showing their best 
model of the reversals in the Milky Way field. 

If we compare figured with subsequent plots, the ongoing magnetic field injections gives an important 
additional mechanism to generate magnetic reversals. In addition to the global scale magnetic field reversals 
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discussed above, this model contains several local reversals which appear to be due to local field injections, 
and which exist for some time. Possibly, such local rev ersals can be compared with reversals suggested to 



We appreciate that the models presented are very oversimplified, but our intention is to show that 
initial conditions can strongly influence the long term evolution of the large-scale magnetic field. With 
"standard" weak initial fields (\B\ « B cq ), any reversals are transient features. It appears important 
that the nonlinearity operates before the large-scale field forms. Of course, confirmation or revision of 
the conclusions of this paper in the framework of 3D mean-field models, or even by direct numerical 
simulations of the microscopic induction equation, performed with realistic hydrodynamical models of 
very early galaxies would be an important development of the theory of galactic dynamos. However this 
obviously is beyond the scope of this paper. 

We are grateful to our colleagues, who work on magnetic field evolution in very young galaxies namely 
T. Arshakian, R. Beck, M. Krause, R. Stepanov, for fruitful discussions which stimulated the writing of 
this the paper. We also thank Anvar Shukurov and an anonymous referee for critical readings of the paper. 
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